Grain boundary pinning and glassy dynamics in stripe phases 
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We study numerically and analytically the coarsening of stripe phases in two spatial dimensions, 
and show that transient configurations do not achieve long ranged orientational order but rather 
evolve into glassy configurations with very slow dynamics. In the absence of thermal fluctuations, 
defects such as grain boundaries become pinned in an effective periodic potential that is induced 
by the underlying periodicity of the stripe pattern itself. Pinning arises without quenched disorder 
from the non-adiabatic coupling between the slowly varying envelope of the order parameter around 
a defect, and its fast variation over the stripe wavelength. The characteristic size of ordered domains 
asymptotes to a finite value Rg ~ Ao e"^/^ exp(|aj/y^), where e ^ 1 is the dimensionless distance 
away from threshold, Ao the stripe wavelength, and a a constant of order unity. Random fluctuations 
allow defect motion to resume until a new characteristic scale is reached, function of the intensity 
of the fluctuations. We flnally discuss the relationship between defect pinning and the coarsening 
laws obtained in the intermediate time regime. 



I. INTRODUCTION 

The motion of topological defects in two dimensional 
smectic phases is studied at a finite distance from thresh- 
old. We focus on the Swift-Hohenberg model of Rayleigh- 
Benard convection and related amplitude equations to 
address the role that non-adiabatic effects play in do- 
main coarsening of a modulated phase, defect pinning, 
and the appearance of glassy behavior. 

Topological defects are often the longest lived modes 
of a non-equilibrium system, with their motion determin- 
ing the longest relaxation times of the structure. Phe- 
nomenological models of defect motion that are based 
on a mesoscopic description have been known for some 
time Such a description, valid for distances much 

larger than the defect core, typically involves time de- 
pendent Ginzburg-Landau equations or their generaliza- 
tions. A few cases have been studied extensively, includ- 
ingdomain coarsening in 0(N) models in nematics 
1^ ^ , and in smectic phases as effectively encountered in 
models of Rayleigh-Benard convection or lamellar phases 
of block copolymers |p|-p^. In the case of modulated 
phases, the motion of a single defect has been widely 
studied within the well known amplitude equation for- 
malism. This method describes the spatio temporal evo- 
lution of the envelope of a base periodic or modulated 
structure [p6|-p0[. The amplitude equation description 
is valid only close to bifurcation points where the spa- 
tial scale of variation of the amplitudes is large or "slow" 
compared with the "fast" period of the base pattern and, 
in the present case, with the extent of the defect core as 
well. 



Far enough from the bifurcation threshold of the mod- 
ulated phase, the separation between slow and fast scales 
no longer holds, and corrections to the amplitude equa- 
tions appear because of the coupling between both scales. 
These corrections are generically referred to as non- 
adiabatic effects. One manifestation of non-adiabaticity 
is that a defect that would be expected to move at con- 
stant velocity from an amplitude equation analysis may 
instead remain immobile or pinned pT|-|2^]. We argue 
below that non-adiabatic effects and defect pinning have 
important consequences for domain coarsening of modu- 
lated phases in two dimensions, and are responsible for 
the formation of glassy configurations. 

Our results complement recent research on glassy prop- 
erties of stripe phases. It has been suggested that systems 
in which long ranged order is frustrated by repulsive in- 
teractions (the latter often leading to the formation of 
stripe phases or other patterns in equilibrium) may in 
fact exhibit the properties of structural glasses. An ex- 
ample are the glassy states recently observed in doped 
semiconductors in a stripe phase Coarse grained 

models with competing interactions of the type used here 
(and also used to study block copolymer melts in lamellar 
phases) have been reintroduced to describe the formation 
of glasses in supercooled liquids Additional equi- 

librium studies of the same models in three dimensions 
based on replica calculations |^ or Monte Carlo simu- 
lations 1^ have been used to argue for the existence of 
an equilibrium glass transition. Structural glasses form 
spontaneously at low temperature without the presence 
of any quenched disorder, and their properties remain in 
general poorly understood. It is noteworthy that coarse 
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grained models exhibiting glassy behavior in the absence 
of disorder are rare, whereas examples of discrete systems 
are known (e.g. Ising models with next-nearest neighbor 
interactions ps] , ^ ). We present here a two dimensional 
study that indicates a dynamical route to the formation 
of glassy configurations in stripe phases. 

We first analyze the motion of a particular type of de- 
fect, namely a grain boundary separating two domains 
of differently oriented stripes. Earlier asymptotic work 
near onset {i.e. in the limit e ^ 0, where e is the dimen- 
sionless distance away from threshold) is extended to the 
region of small but finite e. In Section grain bound- 
aries are shown to move in an effective periodic potential 
of wavelength Ao/2 (where Aq is the periodicity of the 
stripe modulation) and of magnitude that increases very 
quickly with e. Grain boundaries asymptotically pin as 
the driving force for grain boundary motion decreases. 
It is argued that for any finite e an infinite size system 
will not achieve macroscopic long range order dynami- 
cally following a quench. Rather, the characteristic size 
of a domain will not exceed typical value Rg that is pro- 
portional to Xq e~^/^ exp(|a|/-\/e), where a is a constant 
of order unity. 



In Section [II , we incorporate the effect of random fluc- 



tuations and derive the corresponding amplitude equa- 
tions valid for fluctuations of small amplitude. The 
asymptotic motion of a grain boundary can be recast 
as an escape problem in which the effective activation 
barrier is seen to be proportional to the grain boundary 
perimeter. 

Our approach must be considered only qualitative in 
nature because of the scope of the description employed. 
Ginzburg-Landau equations, and more generally ampli- 
tude or order parameter equations (of which the Swift- 
Hohenberg model described below is but one example) 
are only asymptotic, large length scale approximations to 
the physical system they model in the immediate vicinity 
of a bifurcation point. Therefore any short scale phenom- 
ena involved in the description of non-adiabatic correc- 
tions clearly falls beyond their range of validity, at least in 
a systematically quantifiable way. It is nevertheless not 
unreasonable to expect that non-adiabatic effects of the 
sort encountered in order parameter equations will also 
occur in the physical systems which they model. Further- 
more, our results also provide insights into many existing 
numerical studies of these order parameter models, as de- 
scribed below. 



In Section IV, we address the consequences of pinning 



on the domain coarsening that occurs in the intermedi- 
ate time regime following the quench. This subject has 
been the focus of several numerical studies |pl]|-[l5| and, 
more recently, of experimental studies in block copolymer 
thin films |Q and in electro-convection in nematics . 
The results of Sections || and |ll| provide a possible inter- 
pretation of conflicting results in the literature. Previous 
studies of this problem addressed the existence of 

self-similarity during domain coarsening and attempted 



to quantify the time dependence of the linear scale of 
the coarsening structure. The statistical self-similarity 
hypothesis asserts that after a possible transient, con- 
secutive configurations of the coarsening structure are 
geometrically similar in a statistical sense. As a conse- 
quence, any linear scale of the structure (e.g., the average 
size of a domain or grain of like oriented stripes) is ex- 
pected to grow as a power law of time l{t) ^ t^/'^, with z 
a characteristic exponent. Self-similarity is a well known 
feature in systems that order in uniform phases of broken 
symmetry p2| , p3| . However, the determination of z has 
been problematic for stripe phases. Its value appears to 
depend on the quench depth (the value of e), whether 
or not fluctuations are included in the governing equa- 
tions, on the thermal history of the system, and on the 
particular linear scale analyzed. 

Recent work in the limit e showed that coarsen- 
ing is self-similar and that z = 3 |]l5|. The value z — 3 
in that limit can be justified by a dimensional analysis 
of the law of grain boundary motion. We focus here on 
the case of finite e (in practice e > 0.1 for the Swift- 
Hohenberg model), and report a slowing down of phase 
ordering dynamics with increasing e, in agreement with 
the literature. We attribute this behavior to partial pin- 
ning of defects that becomes increasingly important at 
long times as the driving force for coarsening decreases. 
At even longer times, coarsening stops altogether and the 
system reaches a glassy state as the linear scale of the 
structure reaches the critical value Rg{t) computed in 
Section When random fluctuations are incorporated 
in the model, we show that, sufhciently close to onset, 
the value of z remains independent of the intensity of 
the fluctuations, thus verifying the universality implied 
in the self-similarity hypothesis in that region. At larger 
e, we find that fluctuations accelerate ordering kinetics, 
also in agreement with the literature, and that, as ex- 
pected, defect motion is allowed beyond the scale given 
by Rg. At even later times the system orders very slowly, 
possibly logarithmically in time. 



II. NON-ADIABATIC CORRECTIONS AND 
GRAIN BOUNDARY PINNING 

We consider the Swift-Hohcnberg model of Rayleigh- 
Benard convection ||3^ as a prototypical model of a mod- 
ulated phase. The numerical results presented below 
have been obtained from a direct numerical solution of 
the model. The analytic results, on the other hand, follow 
from the corresponding amplitude equation, and hence 
are expected to be of somewhat wider generality. The 
model equation studied here is 
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where ij) is an dimensionless order parameter related to 
the vertical fluid velocity at the mid plane of a Rayleigh- 
Benard convection cell, e is the reduced Rayleigh number 
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{R — Rc)/Rc ^ 1 {Rc is the critical Rayleigh number 
for instability), and ko = 2tt/Xq is the roll wavenumber 
(in Appendix A we outline the connection between this 
model and other coarse-grained models with long range 
repulsive interactions [3^]). 

For < e <C 1, the leading order approximation to the 
stationary solution of Eq. (|l]) is a sinusoidal function of 
wavenumber ko. We focus in this section on a configura- 
tion that contains an isolated grain boundary separating 
two such stationary solutions with mutually perpendic- 
ular wavevectors (Figure [^). The reason for studying 
this perpendicular orientation is the expectation that a 
90° grain boundary is that of lowest energy, and hence 
the prevalent boundary angle in an extended system that 
evolves spontaneously from an initially disordered config- 
uration (see, for example. Figure ^,b). It is known that 
a planar grain boundary separating two regions of uni- 
form ko is stationary [l9|p0[ |. However, we found in ref. 
p6| that a slightly perturbed boundary undergoes a net 
translation with a speed that is a function of the curva- 
ture of the rolls ahead of it. We address in this section 
the extension of the asymptotic results given in that ref- 
erence to small but finite e, and show how corrections 
obtained lead to boundary pinning. 

Near threshold, a 90° grain boundary configuration is 
an approximate solution of Eq. (^ of the form 
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where slow variables are denoted by capital letters and 
are defined as MM 
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(The coordinate x is directed along the normal to the 
reference planar boundary.) 

We recall first some known results for a planar and 
stationary grain boundary in the limit e — > 0, a case 
that was extensively studied in refs. jl^,|2^. The station- 
ary amplitudes {Aq,Bq} are a function only of a;. Aq, 
the amplitude of the rolls parallel to the interface, van- 
ishes as exp{x^/€/Xo) when x — oo, and saturates to 
(4e/3)^/^tanh(a;-v/e/Ao) when x +oo. The behavior of 
the amplitude of the rolls perpendicular to the interface 
is slightly different: Bo{x) - (4e/3)^/2 oc exp(xVe/Co) 
when a; — *■ — cxD and there exists a location x* such that 
Bq(x > X*) to a good approximation. Hence, the 
grain boundary region has a thickness proportional to 
Xo/^/e. It is important to note that at small e the location 
of the grain boundary decouples from the phase of the 
stripes of domain A. Thus, the configuration obtained 
is invariant under any translation of the grain bound- 
ary by a distance Xq (the phase of the stripes remaining 
unchanged) . 



We next derive two coupled equations for the ampli- 
tudes A and B that take into account the possible cou- 
pling between these amplitudes and the phases of the 
stripes. This coupling becomes significant at a finite 
value of e, and hence when there is a large but finite 
separation between the scales {Xa,b,Ya.b} and \x, y}. 
We follow an approach similar to that used in ref. j23| to 
study the motion of a planar front between a hexagonal 
and a uniform phase, or between a hexagonal and a stripe 
phase. The first step is a multiscale analysis, and is stan- 
dard Equation (|l|) is expanded in power series of e, 
as well as the solution ip = e^/^i/'i/2 + eV'i + ^^^'^i^3/2 + ■■■■ 
The leading order solution e^/^Vi/2 is given by Eq. (||). 
At order e"^/^, the solvability conditions for the existence 
of a nontrivial solution for ^3/2 yield the relations that 
A and B must satisfy, 
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with the linear operator C = 1 — dr — k^^ ^(^x^ 
2dydYB + '^(^xOxa)'^- In the limit e — > the functions A 
and B remain constant over one spatial period Aq, and 
therefore the only non vanishing contribution to the in- 
tegrals come from the terms proportional to e*''°^ (resp. 
gjfcoy') within brackets in Eq. (|) (resp. Eq. (|)). This 
standard set of coupled Ginzburg-Landau equations fol- 
lows ||l^,|l6| . It is known, however, that additional non 
perturbative contributions arising from the term '0i/2 ap- 
pear in Eqs. (^) and (^). We focus next on these con- 
tribution and their effect on the relaxation of a slightly 
perturbed grain boundary. 

Integrals of the type J^^^° dx' e™'^''^'^"^^ in Eqs. 
(^ and (^ (where to, n and p are integers) will not 
integrate to zero if the thickness of the grain bound- 
ary profiles along the x direction is finite. (Contribu- 
tions from the direction transverse to the grain bound- 
ary, dy' e'"^''°y'A"BP, wiU be neglected. They are 

typically of the order of B^dyA and hence always smaller 
than the leading analytical terms of the amplitude equa- 
tions.) Terms proportional to e™*^"^ will contribute to 
Eq. (^, and terms proportional to girnkox' +ikoy' 
(^). If we only retain the lowest order term e^/^'0i/2 as 
given by Eq. (||), wc find that only A^ (to — 3) con- 
tributes to Eq. (||), while SA^B (m = 2), as weh as 
3A^B (m — —2, with A the complex conjugate of A), to 
Eq. (|^). Reintroducing the original unsealed variables, 
the generalized amplitude equations read 
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where Fgi, = J dfTgi, is the standard Lyapunov func- 
tional corresponding to the 90° grain boundary. Its vari- 
ational derivatives satisfy [p9 20 
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The last terms in the right hand sides of Eqs. (^) and 
(0) depend on both fast and slow spatial scales, and they 
embody the so-called non-adiabatic coupling between the 
two. Analyzing the effects of these two terms on the 
relaxation of a perturbed grain boundary is the subject 
of remainder of this section. 

We now introduce a small perturbation to the pla- 
nar boundary as shown schematically in Figure |l|. The 
phase of the stripes of domain A is distorted by a uniform 
perturbation of wavenumber g <C /co (and of amplitude 
Sxq <C Aq) in the direction transverse to the stripes. As 
shown in ref. |36j, approximate solutions to Eq. (^)-(0) 
are given by 



B ^ Bo{x- Xgbit)) , 



(10) 
(11) 



where Xgb(t) represents the time dependent position of 
the grain boundary (averaged over y). As already dis- 
cussed in that reference, perturbations to the phase of B 
are of higher order in e. In order to derive a law of motion 
for Xgb it is simpler to neglect the linear relaxation of the 
perturbed rolls, and hence assume that Sxq is constant. 
The amplitude 5xq relaxes exponentially with time but 
the relaxation time of the perturbation is proportional to 
and usually much longer than the characteristic time 
associated with grain boundary motion, Xo/xgb- Further- 
more, as was shown in ref. [ pq] , explicitly considering 
stripe relaxation does not change the law of motion for 
Xgb in any quantitative way. 

Multiply Eq. d) (resp. Eq. (§)) by OtA (resp. dtB), 
add the results and integrate the real part over the system 
area. By using Eqs. (^-(^ and integrating by parts 
the non-adiabatic terms, we obtain the following law of 
motion for the grain boundary 
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where n = Sxqq'^ is proportional to the mean curvature 
of the stripes of domain A, (j) is a, constant phase, and 
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Equation 
in Ref. 
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without the oscillatory term was derived 
in the limit e ^ 0. The coefficient -D(e), 



with dimensions of an inverse length, represents a friction 
term that depends on the static grain boundary profile 
{Aq,Bq\, while the term en"^ in the numerator is pro- 
portional to dy [Tgbix = oo,y) - Tgbix ^ -00, y)]/L, 
where Tgb is the free energy density implicitly defined by 
Eqs. (|)-(|)). The numerator can be understood as the 
leading contribution (in e and k) from an external force 
acting on the grain boundary. This force results from the 
difference in the free energy density J^gb between curved 
stripes on one side, and straight stripes on the other side 
of the boundary. Note the unusual dependence of Xgb on 
a even power of the curvature thus indicating that the 
motion of the grain boundary is such that curved paral- 
lel rolls of higher energy are always replaced by straight 
perpendicular rolls. 

The last term in the right hand side of Eq. (|l^ ) is 
the dominant contribution arising from the non-adiabatic 
terms of Eqs. (|^)-(0). The dimensionless quantity p(e) 
plays the role of the amplitude of a periodic potential 
of period Ao/2 within which the grain boundary moves. 
The major contribution to p(e) comes from the integral 
that contains the term dxBo in Eq. (|lj) since the profile 
Bf)ix) has a steeper variation than ^0(2:) jl^]. Given 
that both amplitudes Aq and Bq are approximately of 
the form y^/(y^a;/Ao), it is easy to show from Eq. (|lj) 
that 
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where \a\ is a constant of order unity, corresponding 
to the pole of the envelopes closest to the real axis in 
the complex plane. Hence, p behaves non-analytically at 
small e, and increases extremely quickly with e. Quali- 
tatively similar results were reported in |^,^ for one 
dimensional fronts between conductive and convective 
states, or between different convective states. 

Equation ( |l2|) shows that for any finite e > a planar 
grain boundary (k — 0) can have only two stationary 
positions per period of the stripe pattern Aq. This effect 
had been observed numerically and reported in ref. [ p6[ , 
with similar findings also given in ref. |po|| . Equation 
( p^ ) also implies that there exists a critical curvature Kg 
below which the grain boundary will remain immobile. 
This critical curvature is given by. 
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where Rg is the associated radius of curvature which di- 
verges non-analyticaUy near onset 

Rg ^ Ao .-V2 exp (^i^) . (17) 

These results have been verified by direct numerical 
solution of the Swift-Hohenberg model with reasonably 
small values of e. The numerical algorithm used has been 
described in |3^Jl^. Briefly, Eq. (^ is discretized on a 
square grid of mesh size Aa; = 1 with 512^ nodes (256^ 
for e = 0.5), and the wavelength is set to Aq = 8Aa;. A 
semi-implicit spectral method is used to iterate in time. 
The initial condition for "0 is a white and Gaussian ran- 
dom flcld with zero average and variance (^^) = e. Typ- 
ical long time configurations which are stationary for all 
practical purposes are shown in Figs. |2| and ||a. These 
figures show the field ip in grey scale. Many topologi- 
cal defects including dislocations, -f 1/2 disclinations and 
several 90° grain boundaries can be identified. Figure 
corresponds to e = 0.5 and two different times t = 10 
and t = 2 X 10^, showing that the order parameter does 
not change beyond t = 10^. Figure ^ corresponds to 
e = 0.4, and the configuration shown remains practically 
constant beyond t = 2.3 x 10*^. 

To further quantify these observations we have com- 
puted the probability distribution function of stripe cur- 
vatures P{K,t). The stripe curvature is defined as k = 
|V • n|, where n is the unit normal to the lines of con- 
stant ip. The curvature k is a slowly varying quantity 
away from defect cores, and only these regions are used to 
compute P(k, t) by the filtering method described in ref. 

Figure ^ shows our results for e = 0.4 and e ~ 0.5. 
In both cases the distribution converges at long times 
towards a limiting curve of finite width, thus indicat- 
ing that asymptotic configurations contain many curved 
stripes and are disordered at large scales, or "glassy". 
This behavior is to be contrasted with that of a coarsen- 
ing system in which P{k, t — > oo) would approach a delta 
function at k = 0. We take P{k — Q,t ^ oo) as a mea- 
sure of the linear scale of the structure or typical domain 
size and compare its value with the pinning radius Rg 
given in Eq. (p^). Figure^ shows the numerical results 
together with Rg multiplied by a (fitted) scale factor ap- 
proximately equal to 4. The pinning radius Rg increases 
extremely quickly with decreasing e, in agreement with 
the numerical calculations for the range of e which we can 
study (computational constraints on system sizes have 
prevented us from investigating the region e < 0.30). We 
have checked that the glassy configurations at long times 
do not result from numerical pinning; the results are not 
modified when the grid spacing is halved to Aa; = Ao/16. 

Although other types of defects (e.g., dislocations and 
-1-1/2 disclinations) may also become pinned, and thus 
contribute to the overall stability of glassy configurations, 
the predominance of grain boundaries over other defects 
seems to be a generic feature of the Swift-Hohenberg 
model (see Figures Ha,b and ref. ) . Furthermore it is 



likely that a similar dependence between the speed of the 
defect and e will hold for the motion of other topological 
defects (except for dislocation climb). Hence we argue 
that a defected configuration of stripes does not macro- 
scopically order following a quench to a finite value of e. 
Asymptotic long time configurations appear to exhibit a 
labyrinthic and partially disordered structure with many 
immobile defects that do not anneal away. These disor- 
dered configurations resemble those of a structural glass 
at zero temperature which lack long range order (trans- 
lational or orientational) order. They become sponta- 
neously trapped in metastable configurations that are 
very different from the configuration of lowest free en- 
ergy (all stripes parallel to each other, or a "crystalline" 
state). 

Not all grain boundaries in a glassy configuration are 
90° boundaries. However, we expect that grain bound- 
aries with a different orientation would be pinned less 
efficiently {i.e. would have a higher value of \a\ in Eq. 
(p^). The reason is that their stationary planar pro- 
file is smoother than that of a 90° grain boundary, and 
therefore non-adiabatic effects are expected to be weaker. 

We finally mention that if both e and kAq are not small 
compared to one, both adiabatic and non-adiabatic terms 
will contain higher order analytic corrections which we 
have not calculated. 



III. MOTION AT FINITE TEMPERATURE 

Given the results of Section ||, it is natural to study the 
effect of random fluctuations added to Eq. (|l|) . Small am- 
plitude fluctuations will allow activated motion of grain 
boundaries, and in general unpinning. We consider in 
this section the stochastic Swift-Hohenberg model 

^ = - + V2)V - V'' + r;(f, t) , (18) 

where ry is a Gaussian and white random noise of zero 
mean and variance 

{ri{r, t)ri{f',t')) = 2F5{r- r ')6{t - t'). (19) 

The noise intensity F is proportional to the (dimension- 
less) temperature according to the fluctuation-dissipation 
theorem. In what follows, F and e are considered as in- 
dependent parameters, although they might be related 
in some particular physical systems. The stochastic 
Swift-Hohenberg model has been used to study hydro- 
dynamic fluctuations near onset of Rayleigh-Benard con- 
vection , and thermal fluctuations of molecular origin 
in lamellar phases of diblock copolymers js^ . 

The stationary states of Eq. ( p^ ) in two spatial di- 
mensions have been studied in refs. ||§,|ll|. Above a 
critical noise intensity F^ (that depends on e), the system 
is disordered (lacks both translational and orientational 
long ranged order). Below Fc a stripe phase with long 
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ranged orientational order but no translational order was 
found. Only at = the system was seen to exhibit 
both translational and orientation long ranged order. In 
what follows we focus on defect dynamics in the range 
< F ^ Fc, so that the local stripe pattern is not very 
distorted. 

We first derive the stochastic amplitude equations for 
a 90° grain boundary. Following Graham I^^Ij 
proximate the effect of the noise on the amplitudes by 
projecting it along the two slow modes of the determinis- 
tic equation and neglecting any contribution arising from 
couplings and resonances between noise and fast variables 
(El|,E2|. We start by writing the random function as, 
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where the slow variables X{Y)a,b are given by Eq. (|^), 
and f/A and f/B are two independent complex random 
processes that satisfy the relations, 
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It is implicit in the decomposition that F is small enough 
so that well defined stripes exist locally. On the other 
hand, F has to be large enough so that f]A and fjB are not 
negligible in the solvability conditions at order e"^/^ . 
Given both assumptions, Equations (|^) and (|^) straight- 
forwardly generalize to 
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We can now estimate the escape rate of a grain bound- 
ary over the potential barrier of Eq. (^2|). In order to 
do so, we need to estimate the projection of the noise 
intensity in Eqs. (|l]) and (H) on the coordinate Xgb{t) 
implicitly defined by Eqs. ( po| ) and ([Tl|). A rough esti- 
mate that is sufficient for our purposes can be obtained by 
using Eqs. ( |lO| ) and ( [Tl| ) as the trial solution of Eqs. ( |2l| ) 
and (|2^). Focusing on Xgb alone ignores possible bound- 
ary broadening because of fluctuations, or roughening. 
Both phenomena will be important for grain boundary 
motion above the pinning point, but their contribution 
is probably less important in the immediate vicinity of 
the pinning transition. By substituting Eqs. ( p^ and 
(|ll|) into Eqs. ^) and (ph, we find 



Xgb 



- 4^ cos(2A:o2:g6 + (/))+ 7^ , (23) 



3fc2D(e) Die) 
with fj a (real) random white Gaussian noise satisfying, 

{fj)=0, (r)(i)77(t')) =2F' 5(<-<'), F' = F/[2D{e)Rgb] 

(24) 

where Rgb is the grain boundary perimeter. As expected, 
the intensity of the fluctuations on the glob al coordinate 
Xgb is proportional to l/Rgb- Equation ( |23| ) is a straight- 
forward generalization of Eq. (|l2|), and is formally anal- 
ogous to the equation that describes the one dimensional 
motion of a Brownian particle in a periodic potential of 
amplitude 2p{e)/[2D{e)ko]. 

Equations (p3[) and (|23) can be recast as 



Xgb 




— cos{2koXgb 



(25) 



The random term ^ is such that (^) = and {£,(t)£,{t')) = 
25(t — t'). We have also used Eq. ( p^ ) to eliminate p{e) 
from Eq. (p3|), and we have defined 



^^0 



2e 



ik^Rg 



(26) 



Consider the situation where grain boundaries are pinned 
at F = 0. Since k < Kg, the first term of the right-hand- 
side of Eq. (^5|) is not dominant and the potential barrier 
that a pinned defect of size Rgb has to overcome is of the 
order of F^/Rg. The stochastic problem is now an escape 
problem over this potential barrier given the intensity of 
the noise term in Eq. (p5|). The Kramers rate of escape 
is given by 



VB- (22) 



r ^ exp — — 



FoR 



gb 



F Ra 



(27) 



Therefore a noise intensity 



F = R 



Ro 



Rgb k^\^ e-l"l/-^ 



(28) 



is required to un-pin a grain boundary of length Rgb . 



IV. SLOW COARSENING DYNAMICS: 
DEPENDENCE ON TEMPERATURE AND 
QUENCH DEPTH. 



We use here the results of Sections || and III to pro- 
vide a possible interpretation of conflicting results con- 
cerning domain coarsening of stripe phases. We recently 
studied this issue by numerically solving the noiseless 
Swift-Hohenberg equation (Eq. (|l|)) in the limit e — > 
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|g5[ . Our numerical results suggested that the character- 
istic scale of the structure (or the typical size of ordered 
domains) increases as i^^^, with z = 3. That value of 
the exponent was interpreted to follow from the domi- 
nant motion of grain boundaries through a background 
of curved stripes. In disordered configurations, the curva- 
ture of stripes is set by a distribution of largely immobile 
-1-1/2 disclinations. According to Eq. (p^, the motion of 
grain boundaries is driven by stripe curvature, and acts 
to reduce the overall curvature by replacing regions of 
curved stripes by straight ones of a different orientation. 
It also reduces the disclination density whenever their 
core region is swept by a moving grain boundary. In the 
limit e <C 1 we computed several measures of the lin- 
ear scale, including moments of P{K,t), moments of the 
structure factor of the order parameter, and the average 
distance between defects. They were all found to become 
proportional to each other, and to grow as a power law 
of time with an exponent 1/3. 

Grain boundary motion as described in Section || was 
used to provide an interpretation for the value z = S. 
Since -fl/2 disclinations generate roughly axisymmet- 
ric patterns of stripes around them, the characteris- 
tic stripe curvature in any given configuration is pro- 
portional to the inverse characteristic distance between 
disclinations. Under the self-similarity hypothesis, the 
distance between disclinations is proportional to l{t), 
hence k ^ l/l{t). If grain boundaries are the class of 
defect the motion of which controls asymptotic coars- 
ening, then the coarsening exponent can be inferred by 
dimensional analysis of Eq. (12). In the limit e ^ 0, the 
oscillatory term in the right hand side of Eq. (|lj) can be 
neglected and we simply have dl/dt cx l^^ or l{t) ~ t-^^^, 
in agreement with the numerical solution of Eq. (|l|) . 

Equation (|l^) shows that this result changes qualita- 
tively further from onset. As e increases the pinning po- 
tential energy barrier p{e) increases extremely fast, and 
important corrections to scaling are to be expected. For 
finite e and short times many defects are present, there- 
fore the characteristic curvature of the stripes is very 
large, and the first term in the right hand side of Eq. 
( p^ ) dominates. As coarsening proceeds, the characteris- 
tic curvature decreases until it reaches the critical value 
Kg given by Eq. (^6|). At that point the typical velocity 
of a grain boundary vanishes, although the system is still 
disordered. Therefore one would expect that coarsening 
would stop when l{t) is of the order of Rg. This is pre- 
cisely the results shown in Fig. ^ with only one adjustable 
parameter (a scale factor relating Rg given by Eq. ( |l^ ) 
to l{t) determined numerically from the distribution of 
stripe curvatures). 

When random fluctuations are considered (F > in 
Eqs. (|l8|)-(P^)), some of the grain boundaries in a frozen 
configuration are expected to resume motion. We argue 
that the structure will continue coarsening until the aver- 
age domain size reaches a new characteristic size Ip > Rg 
that can be estimated as follows. We write a general phe- 
nomenological evolution equation for the domain size l{t) 



directly from Eq. (p5|): 

dl _ 
di ~ 




cos(2fco/ -I- 0) 



(29) 



where we have assumed that, prior to pinning, the various 
length scales remain approximately proportional to each 
other. Recall from Eq. ( pTj ) that F = Fq is required to 
unpin a configuration obtained in the absence of noise, 
for which l{t) Rg. According to Eq. (p9|), coarsening 
proceeds ii F > Fq until a new characteristic pinning size 
is reached given by FqIp/ {FRg) = 1 or 



Ra 



F 



koF- 



.l"l/v^ 



(30) 



After reaching the scale Ip, domains are expected to 
coarsen very slowly by thermal activation. When a grain 
boundary overcomes one pinning barrier, the linear ex- 
tent of the corresponding domain typically increases by 
an amount of order Ao/2. Hence dl/dt ~ Aqt, where r is 
given by Eq. (|2^) with i?gf, replaced by I. Hence domains 
are expected to grow logarithmically in time according to 



l{t) Fln{t/F) for/>Zi 



(31) 



A numerical solution of Eq. ( p^ ) yields results qual- 
itatively consistent with those presented above. Figure 
^ shows a configuration of the order parameter field ip 
obtained for F = and e = 0.4 starting from random 
initial conditions. The configuration shown corresponds 
to very late times t = 2.3 x 10^ at which point all de- 
fects are practically immobile, and domain growth has 
stopped. We then set F = 0.00318, and the integration 
is continued. The order parameter configuration t = 10^ 
time units later is shown in Fig. |^b. The average domain 
size has increased substantially. Many grain boundaries 
have a 90° orientation (like in Fig. ||a), and roughening is 
limited or nonexistent. We have determined the average 
domain size I from the probability distribution function 
of the quantity ( = t/;'^ + iyipY /k^. Figure ^ shows 
the probability distribution function corresponding to a 
perfectly ordered configuration, as well as to partially 
disordered configurations. The inverse linear scale 1/Z, 
proportional to the defect density pd, is extracted from 
the difference between these curves, as detailed in Ap- 
pendix ^ As shown in Fig. |^ domain growth is very 
slow, possibly logarithmic, although a precise check of 
this behavior is problematic. 

Figure || displays the evolution of the defect density 
Pd{i) as a function of time, starting from random ini- 
tial configurations. For reference we also show the case 
F — Q. Increasing the value of e leads to smaller effective 
exponents, whereas increasing F has the opposite effect. 
For sufficiently small e, we find z — Z independent of 
the value of F. The two bottom curves correspond to 
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systems that are close enough to onset, and hence either 
Rg or Ip is very large compared with the linear size of 
the system. We show our results for e = 0.04 (averaged 
over 40 independent runs) and for e = 0.15 (averaged 
over 15 independent runs). The solid line closest to these 
two curves has a slope of —1/3. The downward deviation 
from linearity at long times at e = 0.04 is a typical mani- 
festation of finite size effects (this long time behavior and 
its dependence on the system size was studied in detail 
in ref. ||). 

With increasing e and/or decreasing F, pinning be- 
comes more pronounced as evidenced the lower effective 
slopes of the three upper curves in Fig. || (the results 
are averages over six independent runs, each curve corre- 
sponding to the same value e = 0.4). The top curve cor- 
responds to a system without fiuctuations for which pd 
was computed with the method described in ref. fl^ j . The 
density starts decaying roughly as an inverse power law, 
with an effective exponent much smaller than —1/3 (the 
top solid fine has a slope of —1/5), and after a crossover 
saturates at long times indicating pinning. When small 
amplitude noise is added (curve below denoted by dia- 
monds), the initial behavior is similar to that of = 0, 
and the decay rate also slows down considerably at long 
times (where we would predict logarithmic growth, i.e. 
Pd l/lii(0)- The curve below, denoted by plus signs, 
corresponds to a noise intensity three times larger than 
the previous case. Its initial decay is slightly faster (it 
can be fitted with an effective exponent 1/ze// — —0.23 
as shown with the solid line in the figure), and the up- 
wards deviations at long times are less pronounced. This 
behavior is in qualitative agreement with the the expecta- 
tion that defects overcome pinning barriers more readily 
at higher noise intensities and pinning is postponed to 
longer times when Z ~ Zp. However, the effective ini- 
tial decay is slower than t^^/"^ which we interpret as a 
crossover effect resulting from non-adiabaticity. 

In summary, our results for large e are in agreement 
with earlier numerical results performed at e = 0.25, 
showing that coarsening laws are very slow and depend 
on the presence of thermal fluctuations. We argue here 
that a coarsening exponent can be properly determined 
only in the small e limit, where the phase ordering kinet- 
ics is self-similar. Our results support that the exponent 
2; = 3 is independent of F for sufficiently small e, when 
pinning effects are negligible {Rg much larger than the 
linear size of the system). 

V. CONCLUSIONS 

We have shown that the Swift-Hohenberg model of 
Rayleigh-Benard convection exhibits glassy properties in 
spatially extended systems. In the absence of fluctua- 
tions, and following a parameter quench across thresh- 
old, random initial configurations do not evolve into com- 
pletely ordered states, a single plane-wave or crystalline 



state. Instead, they reach disordered metastable con- 
figurations in which topological defects, mainly grain 
boundaries and disclinations, fail to annihilate and re- 
main with finite density. It appears that the formation 
of these glassy configurations in a quenched disorder- free 
system can be accounted for by the finite separation be- 
tween "fast" length scales of the structure (associated 
with stripe periodicity), and "slow" scales (associated 
with the extent of defect envelopes) . Since at a finite dis- 
tance from threshold the ratio between these two scales 
is finite, non-adiabatic effects lead systematically to de- 
fect pinning in an infinite system. Fluctuations allow 
unpinning and a certain amount of "crystallization" , al- 
beit through an asymptotically slow activated motion of 
grain boundaries and other defects. 

The present framework is far too simple to be used in 
the prediction of a glass transition temperature, if such a 
transition exists. In some respects, the situation just de- 
scribed is instead very similar to that of domain growth 
in random fields in dimension larger than two p^ , p5[ . 
There, domain walls separating magnetized domains are 
pinned by fixed impurities and evolve by thermal acti- 
vation to other more favorable configurations. The phe- 
nomenological pinning energy of a domain of size R grows 
as Ti?^, were 6 depends on the problem considered, yield- 
ing an escape rate given by r ~ exp(— T^R^/fc^r), equa- 
tion that is formally analogous to Eq. ( p7| ) with 9 — 1. 
Two crucial differences are that our system is glassy even 
in two dimensions, and that defects do not need any dis- 
order to become pinned. 

The consequences of defect pinning on the interme- 
diate time regime corresponding to domain coarsening 
have also been investigated. A universal coarsening ex- 
ponent can be determined close to threshold only, where 
we obtain z — Coarsening stops when the linear 
size of the system is larger than the characteristic do- 
main size for pinning. In this situation, an intermediate 
crossover regime is anticipated with lower effective coars- 
ening exponents, as is observed in numerical solutions of 
the model. Crossover effects induced by pinning can be 
reduced by either increasing the intensity of the fluctua- 
tions or approaching threshold. 

We note that some of our conclusions as well as our 
interpretation of the numerical results are based on the 
analysis of a particular type of defect, namely a grain 
boundary separating two domains with differently ori- 
ented stripes. We think it is likely that similar non- 
adiabatic corrections to defect motion will appear for dis- 
location glide or disclination motion, leading to similar 
non perturbative corrections in e to the speed of the de- 
fect. 

We believe more generally that pinning through non- 
adiabatic effects is likely to be a feature of a wide variety 
of pattern forming systems, and is not limited to the par- 
ticular model treated here. Block copolymer melts, for 
instance, provide an interesting case in which the results 
obtained could have practical implications (see Appendix 
A for a summary of the relevant equations and their re- 
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lationship with the model studied here). We also men- 
tion here that results qualitatively similar to ours have 
been reported for a model with competing interactions 
(describing ferromagnetic films), that is defined by the 
equations of Appendix A with a different form of the 
Green's function G There, frozen polycrystalline 

configurations of stripe patterns were observed for deep 
quenches as well, whereas the system could reach an or- 
dered state for shallow quenches. This same model was 
also able to predict the formation of a frozen phase com- 
posed of poh[disperse droplets with a near-hexagonal ar- 
rangement p7[ | , as previously observed in experiments on 
a Langmuir monolayer p8| . However, the pinning mech- 
anism involved in this last case is probably different than 
the one discussed in the present paper since the patterns 
are no longer locally periodic. Nevertheless, we would ex- 
pect that our main conclusions can be readily extended 
to other systems with periodic structures such as hexag- 
onal patterns [|^,|9i. 
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APPENDIX A: MEAN-FIELD MODEL OF A 
SYMMETRIC BLOCK COPOLYMER MELT 

We briefly recall in this appendix known results about 
the relationship between the mean field description of a 
block copolymer melt, and the amplitude equation for 
Swift-Hohenberg model (Eq. (|l])) at first order in e. The 
dynamics of micro-phase separation of block copolymers 
is often modeled by a time-dependent Ginzburg-Landau 
equation for a conserved order parameter ]35|,|50|, 



gV^(r,t) 
dt 



6F 



(Al) 



where 



F ^ I dr 

B 
'~2 



drdr'i^{f,t)G{f,f')^{r',t). (A2) 



G is the Green's function of the Laplacian operator 
V^G(r, r ') — —6{f— r') and M a constant mobility 
or Onsager coefficient. The scalar order parameter ?/; 
is the local monomer concentration difference between 
the two chemical species. Following Ref. we set 
r = 2 + e,M = f,Ar = f/fcg and B = fcg. Two inde- 
pendent parameters fco and (small) e remain. Equation 
(Al) reduces to 



(A3) 



where ipoc is the boundary condition at infinity. In most 
studies, it is customary to set ipoo = {^), the spatial av- 
erage of ip over the sample. We introduce the amplitude 
A of slightly modulated waves through, 



I 



(A4) 



A multi-scale analysis of Eq. (A3) in the limit e <C 1 
was conducted by Shiwa ||5l|. Setting M = I /kg, the 
resulting equation for the amplitude is 



dA _ 4 
5^ "'^+^2 



A--\A\'A, (A5) 



which is identical to the amplitude equation of the Swift- 
Hohenberg model. Note that the only effect of the con- 
servation law on the local part of the free energy (the 
Lap lacian operator in front of the square bracket in Eq. 
( |A3| )) is a renormalization of the mobility M. The quan- 
tities e and kg defined above play the same role as the 
same coefficients in the Swift-Hohenberg model {i.e. the 
dimensionless distance to threshold and the dominant 
wavenumber of the structure, respectively.) 



APPENDIX B: CALCULATION OF THE DEFECT 
DENSITY IN THE PRESENCE OF 
FLUCTUATIONS 

Computation of the domain size from the probability 
distribution of stripe curvature is delicate in the presence 
of noise. We have used a different method than that used 
for F = 0. We introduce an effective squared amplitude 
C{r, t) by 



C = V-' + (v^)Vfco • 



(BI) 



For a perfectly ordered system consisting of a plane- wave 
solution of the Swift-Hohenberg equation and F = 0, 
the probability distribution function of C is a delta func- 
tion at Coo = 4e/3. When F > the probability dis- 
tribution function of Q even for a plane wave poo (C) is 
broader because of "phonon" excitations. The function 
Poo'' is plotted in Figure ^ (with dotted lines), and is 
then used as a reference curve for a fixed F. In a par- 
tially disordered configuration, the presence of defects 
and curved stripes further broadens the probability dis- 
tribution function to p^^\C,,t) (solid lines of Figure H). 
The difference between the two curves is related to the 
degree or disorder beyond small fluctuations away from 
a perfectly ordered structure. We define the defect den- 
sity by pd{t) = Max{p(^)(C), C} - Max{p(^) (C, t), C}- 
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Since grain boundaries are seen to be the major contri- 
bution to defect density, one can introduce a characteris- 
tic length scale p'^^, which is further identified with the 
characteristic size, or domain size I. In Figure |^, I has 
been normaUzed so that l{t = 0) = Rg (at the beginning 
of the heating process) where Rg is computed from the 
probabiUty distribution function of stripe curvatures at 
F = 0. 



[1] K. Kawasaki, Prog. Theor. Phys. Suppl. 79, 161 (1984). 
[2] H. Brand and K. Kawasaki, J. Phys. A 17, L905 (1984). 
[3] G. Mazenko, in Formation and interaction of topologi- 
cal defects, edited by A.-C. Davis and R. Brandenberger 

(Plenum, New York, 1995). 
[4] A. Bray, in Formation and interaction of topological 

defects, edited by A.-C. Davis and R. Brandenberger 

(Plenum, New York, 1995). 
[5] I. Chuang, B. Yurke, A.N. Pargellis, and N. Turok, Phys. 

Rev. E 47, 3343 (1993). 
[6] H. Toyoki, Phys. Rev. E 47, 2558 (1993). 
[7] M. Zapotocky, P.M. Goldbart, and N. Goldenfeld, Phys. 

Rev. E 51, 1216 (1995). 
[8] N. Mason, A.N. Pargelis, and B. Yurke, Phys. Rev. Lett. 

70, 190 (1993). 
[9] Y. Oono and M. Bahiana, Phys. Rev. Lett. 61, 1109 

(1988). 

[10] K.R. Elder, J. Viiials, and M. Grant, Phys. Rev. Lett. 

68, 3024 (1992). 
[11] K.R. Elder, J. Vinals, and M. Grant, Phys. Rev. A 46, 

7618 (1992). 

[12] M.C. Cross and D.I. Meiron, Phys. Rev. Lett. 75, 2152 
(1995). 

[13] Q. Hon, S. Sasa, and N. Goldenfeld, Physica A 239, 219 

(1997) . 

[14] J.J. Christensen and A.J. Bray, Phys. Rev. E 58, 5364 

(1998) . 

[15] D. Boyer and J. Vinals, Phys. Rev. E 64, 050101 (2001). 

[16] P. Manneville, Dissipative Structures and Weak Turbu- 
lence (Academic, New York, 1990). 

[17] M. Cross and P. Hohenberg, Rev. Mod. Phys. 65, 851 
(1993). 

[18] E.D. Siggia and A. Zippehus, Phys. Rev. A 24, 1036 
(1981). 

[19] P. Manneville and Y. Pomeau, Phil. Mag. A 48, 607 
(1983). 

[20] G. Tesauro and M. Cross, Phil. Mag. A 56, 703 (1987). 

[21] Y. Pomeau, Physica D 23, 3 (1986). 

[22] D. Bensimon, B. Shraiman, and V. Croquette, Phys. Rev. 

A 38, 5461 (1988). 
[23] B.A. Malomed, A.A. Nepomnyashchy, and M.I. Tribeh 

sky Phys. Rev. A 42, 7244 (1990). 
[24] M.-H. Julien et al, Phys. Rev. Lett. 83, 604 (1999). 
[25] D. Kivelson et al, Physica A 219, 27 (1995). 
[26] J. Schmalian and P.G. Wolynes, Phys. Rev. Lett. 85, 836 

(2000). 



[27] M. Grousson, G. Tarjus, and P. Viot, Phys. Rev. Lett. 

86, 3455 (2001). 
[28] J.D. Shore, M. Holzer, and J.P. Sethna, Phys. Rev. B 46, 

11376 (1992). 

[29] M.E.J. Newman and C. Moore, Phys. Rev. E 60, 5068 
(1999). 

[30] C. Harrison et al, Science 290, 1558 (2000). 
[31] L. Purvis and M. Dennin, Phys. Rev. Lett. 86, 5898 
(2001). 

[32] J. Gunton, M. San Miguel, and P. Sahni, in Kinetics 
of first order phase transitions. Vol. 8 of Phase Transi- 
tions and Critical Phenomena, edited by C. Domb and 
J. Lebowitz (Academic, London, 1983). 

[33] A. Bray Adv. Phys. 43, 357 (1994). 

[34] J. Swift and P.G. Hohenberg, Phys. Rev. A 15, 319 
(1977). 

[35] L. Leibler, Macromolecules 13, 1602 (1980). 

[36] D. Boyer and J. Viiials, Phys. Rev. E 63, 061704 (2001). 

[37] G. Ahlers, M. Cross, P. Hohenberg, and S.Safran, J. 

Fluid Mech. 110, 297 (1981). 
[38] P.C. Hohenberg and J.B. Swift, Phys. Rev. E 52, 1828 

(1995). 

[39] J. Toner and D.R. Nelson, Phys. Rev. B 23, 316 (1981). 
[40] R. Graham, Phys. Rev. A 10, 1762 (1974). 
[41] F. Drolet and J. Vinals, Phys. Rev. E 57, 5036 (1998). 
[42] F. Drolet and J. Vinals, Phys. Rev. E 64, 026120 (2001). 
[43] F can be estimated to be of the order F = iie^''*, with 

?i < 1. See also Ref [0. 
[44] T. Natterman and J. Villain, Phase Transitions 11, 5 

(1988). 

[45] J.-P. Bouchaud, L. Cugliandolo, J. Kurchan, and M. 
Mezard, in Spin Glasses and Random Fields (World Sci- 
entific, Singapore, 1998), p. 161. 

[46] C. Sagui and R.C. Desai, Phys. Rev. E 49, 2225 (1994). 

[47] C. Sagui and R.C. Desai, Phys. Rev. E 52, 2822 (1995). 

[48] M. Seul, Europhys. Lett. 28, 557 (1994). 

[49] K.R. Elder, M. Ka takowski, M. Haataja, and M. Grant, 
|cond-mat/0107384 

[50] T. Ohta and K. Kawasaki, Macromolecules 19, 2621 
(1986). 

[51] Y. Shiwa, Physics Letters A 228, 279 (1997). 



10 




FIG. 1. Schematic grain boundary configuration separating two domains of stripes A and B of the same periodicity 
(|fco| = l^o'l — '^o)- The stripes of domain A are weakly curved by a transverse modulation of wavenumber q <C fco. Sxq 
represents the magnitude of the phase modulation. 




FIG. 2. Glassy configurations obtained by numerical solution of the Swift-Hohenberg model with random initial conditions. 
Dimensionless times shown are (a), t — 10000; and (b), t — 20000. Here e = 0.5 and the system has 256^ grid nodes. 
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FIG. 3. Probability distribution function of stripe curvatures, P{K,t), after a quench at e = 0.5 (dotted lines) and e = 0.4 
(solid lines), averaged over 10 and 6 independent runs respectively. For e = 0.5 the figure shows the curves obtained at times 
t = 10*, 5 X 10* and 10^ and for e = 0.4 at times t = 6x 10*, 1.2 x 10^ and 2.3 x 10^ 
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FIG. 4. Characteristic asymptotic grain size following a quench as a function of e. (a), estimate given by Eq. (^^. (b), 
numerical value of P{k = 0,t = oo) (symbols) compared also with Eq.(|l^) multiplied by one fitted scale factor (solid line). 
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FIG. 5. (a) Near stationary configuration obtained after a quench at e = 0.4 and in the absence of fluctuations F = (the 
time shown is t = 2.3 x 10^, and the system size includes 512^ nodes), (b) New structure obtained after taking the configuration 
shown in (a) as an initial condition and further integrating the model equations with F — 0.00318 for a period of 10''' time 
units. At this time, any boundary motion is very slow. 
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FIG. 6. Probability distribution function of ( defined by Eq. (Bl) with ip the solution of Eq. (jl^), for F — 0.00636 and 
e = 0.4. The dotted line corresponds to a single plane wave with superimposed fluctuations, while the solid line correspond to 
disordered configurations obtained from random initial conditions (at times t = 5 x 10'^, 10'* and 10^, respectively). 
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FIG. 7. Characteristic domain size as a function of time at e = 0.4 and F 7^ as indicated in the figure. The initial condition 
at time t = is a glassy configuration obtained from a previous run with F = 0. 
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FIG. 8. Defect density in arbitrary units as a function of time for several values of e and F. The straight solid lines are 
guides to the eye with slopes, from bottom to top, —0.33, —0.33, —0.23 and —0.20. 
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